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We illustrate the general scheme of the Sum Rule (SR) method using 
2D Quantum Harmonic Oscillator (2DQHO) as a toy model. We introduce 
correlator, related to Green function of 2DQHO, and describe the property 
of Asymptotic Freedom for 2DQHO. We explain how the duality conception 
allows one to describe excited states. Finally we present numerical results 
and extract some lessons to learn from our exposition. Then we switch to 
the QCD and show that QCD SRs supply us the method to study hadrons 
in non-perturbative QCD. Here the main emphasis is put on the pion and 
its distribution amplitude and form factors. 

1. Quantum-mechanical toy model 

Two-dimensional oscillator with potential V(r) = mu 2 r 2 /2 is the sim- 
plest system with confinement. We select this particular case D = 2 because 
all formulas greatly simplify, for example, energy levels and wave function 
values in the origin are 

E n = (2n + l)u>; |^„(0)| 2 = — . (1.1) 

7T 

We will consider the regular quasi-perturbative method of Sum Rules (SR) 
to determine energy Eq and | -00 (0) 1 2 of the ground state. We follow here 
partially the lectures pQ. 

The general scheme of the SR method j2] can be easily understood on 
the example of a correlator M(p>), which has the spectral expansion: 

M spec (//) = |^o(0)| 2 e~ Eo ^ + "higher states" . (1.2) 
Suppose that we can construct the perturbative expansion of this correlator: 

2n 

M pert (/i) = Mo(/x) + Y, C 2n^ , (1-3) 
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where Mo(m) corresponds to "free movement" and has the spectral repre- 
sentation: 

M (aO= / Po(E)e- E '»dE. (1.4) 
J o 

The sum rule - it is simply 

M spec (p) = M pcrt (p). 

Usually it appears that higher state contributions can be well approximated 
by "free states" outside interval (0, So). As a result we have SR in the form: 



fSo 2 .4 

|^o(0)| 2 e-^= / p (s)e- s ^ds + C 2 ^ + C 4 ^ + ... (1.5) 
Jo P P 



rS w 2 4 

p (s)e- s ^ds + C 2 — + C A - 
l yr P 4 

Our aim: to determine [^(O)) 2 and Eq from this SR by calculating spectral 
density po (E) of "free particle" and coefficients C 2n by demanding stability 
of this SR in variable p G [p~L,pu] with appropriate value of So. 

In order to select M(p) with these properties for our 2DQHO let us 
consider 2-time Green function 1 

G(0, 0|x, t) = J2 r k (x)M0)e- iEkt (1-6) 

fc>0 

which is the probability amplitude for the transition (x = 0, t = 0) — > (x, t). 
To get M(p) we put x = 0, t = 

M(p) =G(0,0|0,1AX) = J2\ip k (0)\ 2 e - Ek ^ = M s P ec ( y p). (1.7) 

fc>0 

In our case (^(O)] 2 = rnuj/ir, so we have the exact result for our M(p): 

^ 2ir sinh (uj/ ' p) ^ 

How fast is convergence of spectral expansion (|1.7|) for M(/j)? Let us es- 
timate it for p = uj. Exact result ()1.8|) gives us M(cj) = (muj/2-n) • 0.851, 
whereas numerically 

777- UJ 

M spec (uj) = — (0.736 + 0.100 + 0.013 + 0.002 + ...). 
2tt 



Due to the fact that i/>fc(0) = for all states with L > 1 only S-states contribute to 
this sum. 
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Ground state contributes 86%, first excitation - 12%, while the second - 
only 1.5%. 

Now we turn to the perturbative expansion of M(/x) in powers of (o>//z): 



M pert (fi) 



1 



+ 



7 UJ 



31 uj 



6 



+ ... 



2vr V' 6// 2 360 ^ 15120 ju 6 
Here m///2vT corresponds to the Green function of free particle: 



,1/ 



Numerically at [i = uj we have 



muj 



free 



mil 



(1.9) 



(1.10) 



M pert (u;) = — (1 - 0.167 + 0.019 - 0.002 + . . .) . 

First correction specifies free result by 17%, while the second - by 3%. This 
perturbative expansion can be rewritten 



M(/i) - Mp{n) 
M (/i) 



' 6/j, 2 + 360 /i 4 



31 u b 
15120// 



+ 



(1.11) 



that means Asymptotic Freedom: M(/x) behaves like Mo(fi) at large /i>w! 
And it is violated by power corrections of the type (w 2 //i 2 ) . 

We see in Fig.^that for large [i asymptotic freedom works well: M(jx) ~ 
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Fig. 1. Left panel: Exact correlator M{ji), l|1.8|l . is shown by the red solid line, ground 
state contribution only, Mg pcc (fj,), - by the green dashed line, and Mq([i) + 0(u 2 / jj, 2 ) 
- by the blue dashed line. Right panel: Ground + 1-st excited state contribution only, 
^o+T(m); i s shown by the green dashed line, whereas Mo(/x) + 0(w 4 //i 4 ) - by the blue 
dashed line. 

Mo(/u), but we need more and more resonances to saturate M(/x). For 
small ji in spectral part survives only ground state \tpo\ 2 e~ E °^ ,J ', but the 
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perturbation expansion breaks down. In order to model higher resonances 
let us consider the spectral representation of our correlator M(p): 



M spec (/i) = — e~ Ek/l " = rp° SC {E) e~ E/ ^ dE . 
k>o n ^° 



Here the spectral density is just the sum of ^-functions: 

muj 

7T 



P 



\E) = Y J — KE-E k) . 

L * IT 



k>0 



Analogously we have integral representation for free correlator: 



(1.12) 



(1.13) 



(1.14) 



rn 



where po(E) = — . Asymptotic freedom dictates global duality of these two 
densities (term "global" is related with integration over the whole spectrum) 

/ p osc (E)dE= / Po (E)dE. (1.15) 
Jo Jo 

At first glance these spectral densities have completely different behaviour, 
see Fig. |2 But we have very interesting relations between 2fcu;-partial inte- 



p(E)k 



p(E)k S(E - E ) 



Po(E) 



* — * 



Po(E) 



2u) 3w 



7lo 



En Srt 



Fig. 2. Left panel: The exact spectral density p osc (E), is shown by red solid 

vertical lines, imitating S- functions, whereas free spectral density, po(E), - by the blue 
solid line. Right panel: The phenomenological model p Tnod (E), 11.181 . is shown by red 
solid vertical line, imitating 5-function, and by the blue solid line, corresponding to po(E) 
and starting from the threshold So. 

gral moments of this dual densities, namely, {E N ) 2 koj = J2kuj +2U; E N p(E) dE: 



,2(fc+l)u, 

p osc (E)dE = — 



2(k+l)w 



2kui 
2(Jfc+l)w 



Po(E)dE; (1.16) 



Ep {E)dE. (1.17) 



2ku 
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For N > 2 we have approximate relation (E N )° S ^ = {E N )l kuJ [l + O (N 2 /k 2 )] . 
So, we have duality between each excited resonance in oscillator and free 
particle in some spectral domain. That means local duality. 

Now we can model higher state contributions by "higher states" = "free 
states" outside interval (0, So) or: 

p™*(E) = |V> (0)| 2 5 (E - E ) + po(E) 9 (E - S ) (1.18) 

and this gives us 

/•oo 

M mod (p) = |^o(0)| 2 e~ E °^ + / po(s) e-W" alE . (1.19) 

JSn 



After all we have the following SR: 

rS 



\ip (0)\ 2 e" Eo ^ = / p (E) e~ E ^alE + power corrections, (1.20) 
Jo 

or, equivalently, with ^o(O) = V'o(O) \J-k/uj: 



2 



, R()(0) |2 e -^ = £ |l - e"^ - ^ + . . . j . (1.21) 

-l. 



We also have a daughter SR produced from (|1.21|) by applying d/dfj, 

|* (0)| 2 £ e-^ = ^ j 1 " (i + 7) + ^2 + • • j • ( L22 ) 



If we divide (|1,22|) by (|1.21j) we obtain SR for E{ 







p _ , l-(l + 5 /^e-^ + a; 2 /(6^) + ... 

l-e-^- w V(6^) + ... ■ (L23) 

The strategy of processing these SRs is: 

• To determine Eq ~ Eo(So,fj,) by minimal sensitivity to variation of 

€ [/^l; //tj] at appropriate So; 

• To determine | x I , o(0)| 2 ~ *S?q(So,Eq, (j,) by minimal sensitivity to vari- 
ation of fi at appropriate So- 

How we should determine the fidelity window [^l;Mu]? Power corrections 
are of the type (ojj 1 [i) 2n and they are huge at \i <C u. We demand: 

\cM = V ^rrv~ ^ °- 33 for a11 ^ > ml • (1.24) 

Higher states are not suppressed by e~ Ek '^ ~ 1 at large /i>w. We demand: 
A H .s. (ji) = I" e~ E ^dE < 0.33 for all M < Mu . (1.25) 



G 
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Fig. 3. Left panel: Determination of the fidelity window using criteria 11.24H and 11.251 1. 
Right panel: We show here the 1. h. s. of Eq. (11.211 with Eo fixed at the exact value 
Eo = uj as a function of /x (red solid line). Black dashed line corresponds to the position 
of the exact value of |*l/o(0)| 2 . Positions of /xl and /xu are shown by green vertical arrows 
in both panels and fi is displayed in units of uj. 

Then the fidelity window is /jl < /i < f^v'- Only for /i inside it is reasonable 
to demand minimal sensitivity of SRs to variations in /j! 

Let us first consider SR setup with fixed Eq: We fix the energy of the 
ground state to the exact value, Eq = uj, and obtain the following fidelity 
window, /iL = 0.73 a; and /iu = 1.80 a;, see Fig. 01 left panel. The result 
| (0)| 2 = 0.99 is obtained with only 2 power corrections included and for 
So = 2.08a; (the exact value is | x I , o(0)| 2 = 1), see Fig. £3 right panel. 

Now we consider SR in the complete setup, that means that we determine 
the energy of the ground state from the daughter SR (jl.23f) . We take into 
account 3 power corrections and obtain the fidelity window [0.74 a;; 1.8 a;] 
and E = 0.98a; for So = 1.88a;, see Fig. HJ right panel. 

Our conclusions about SRs in quantum mechanics can be summarized 
as follows: 

• SRs give i?o and |V>o(0)| 2 with accuracy not worse than 10% ; 

• Main source of the error - crude model for the spectral density of 
higher states: even taking into account 10 power corrections we obtain 
Eq = 0.95a;, S = 1.79a;, and |^o(0)| 2 = 0.89; 

• But: If we know Eo = 1 exactly (say, from Particle Data Group), 
then accuracy can be twice higher: with taking into account 2 power 
corrections we obtain So = 2.08a; and ^(O)! 2 = 0.99 ! 

• In QCD spectral density more close to perturbative! 
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2. QCD SRs: Way to study hadrons in nonperturbative QCD 

In QCD we have a big problem: nobody knows how to analyze bound 
states. The method of QCD SRs allows us to calculate properties of hadrons 
(masses, decay constants, magnetic moments) without considering hadroniza- 
tion or confinement issues. It was invented in 1977 by Shifman, Vainshtein 
& Zakharov (ITEP) [3] in order to describe J/^-meson, the cc-system, dis- 
covered in 1974 in e + e _ -annihilation at SPEAR (SLAC) and, in parallel, in 
p + i?e-collisions at BNL. In 1979 this method was applied to describe light 
hadrons in massless QCD j2]. 

Main idea: to calculate correlators of hadron currents (0|T [Ji(x) ^(O)] |0) 
by two approaches and to obtain the SR as the result of the matching. We 
start with the Fourier transformed correlator of two hadron currents 

n (Q 2 ) Lorentz 12 = i J e iqx [(0|T [J x {x) J 2 (0)] |0>] dx , (2.1) 

where Lorentzi2 includes all Lorentz structures, so that IT (Q 2 ) is scalar, 
and use the dispersion integral representation 



U(Q 2 ) = I Pl2 + "subtractions" . (2.2) 

s + Q z 



o 

Then we apply Borel transform defined as 
$(M 2 ) = B{Q 2 -» M 2 )II(Q 2 ) = lim 1 " ' 



Tin) 



in 



dQ 



(2.3) 



1.1 


U,(5 




S = 1-88 w 
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Fig. 4. Left panel: We show here the 1. h. s. of Eq. 1)1.2^ . Eo(So, fJ-)/u, as a function 
of fi in units of ui (red solid line). Black dashed line corresponds to the position of exact 
value of Eo- Right panel: The 1. h. s. of Eq. II. 2111 is shown as a function of /i in 
units of u (red solid line) . Black dashed line corresponds to the position of exact value of 
l*o(0)| 2 . 
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Here we list the most important examples: 



n(Q 2 ) 


C = CONST C log (Q 2 /fi 2 ) 


1/Q 2n 1/ (s + Q 2 ) 


$(M 2 ) 


-c 


l/(T(n)M 2n ) e~ s/M2 /M 2 



For us the most important are the first and the last columns in this table: we 
see that Borel transform kills "subtractions" and suppresses "higher states" 
(by the factor exp(— s/M 2 ) in the integrand) in (12,21) : 



B q ^m> [U(Q 2 )} = * (M 2 ) = J Pl2 (s) e~ slM2 ^- 2 (2.4) 

o 

In the 1st approach we apply the Operator Product Expansion (OPE) with 
account for quark and gluon condensates in QCD vacuum to obtain 

if, (r>2\ a (rfl\ _±_ {(asMGG) a s (qq) 2 , 

* (Q ) = ^pert (Q ) + C GG + Cgg M g • (2.5) 

Here ((a s /ir) G a ^ u G a ^ u ) = 0.012 GeV 4 , a s (qq} 2 = 0.0018 GeV 6 has been 
determined in j2j and till nowadays are practically the same [3]. 
The 2-nd approach uses phenomenological saturation of spectral density by 
hadronic states, see (j3.16j) . 

Phad (s) = f%5 (s - m 2 h ) + pp ert (s) 6(s- s ) . (2.6) 

Our model is the ground state h + continuum, which starts from threshold 
s = sq. Then the SR is 

so 

Jh e h - P P ert(s)e as + cqc J^2 h c m M 4 ■ y 2 - 7 ) 



Our aim: to determine /| and from this SR by calculating spectral 
density p peT t(s) and coefficients cqg and Cq q and by demanding stability of 
this SR in variable M 2 6 [M-j^My] with appropriate value of Sq. 

3. Mesonic currents in QCD 

Let us now write down the currents related to ^-mesons in QCD: 

AV: J^(x) = u(x)j^ 5 d(x) ; J^ 5 (x) = d(x)^ ^u(x) ; (3.1) 
PS: Jb(x) = iu(x)^d(x) ; j\ (x) = i d(x)^u(x) . (3-2) 
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Note that Dirac equation iDq(x) = m q q(x) gives us the relation: 

d^J^x) = (m u + m d ) J 5 (x) . (3.3) 

The decay constant f n of the physical pion 7r(P) is defined via 

(0\j^(0)\7r(P))=iUP fM . (3.4) 

It was measured in the decay tt — ► pv^ to be / w = 132 MeV. Eq. ()3.3|) then 
gives 

(0 J 5 (0) vr(P)) = J?r » , (3.5) 
m u + m d 

meaning that the pion reveals itself both in the axial and pseudoscalar 
currents! 

Currents related to /o^-mesons in QCD are 

Jp{x) = u(x)j^d(x) ; J ] p {x) = d(x)^ p u(x) . (3.6) 

The decay constant f p of physical p^(P, e)-meson with polarization e and 
momentum P, satisfying (Pe) = and (e,e) = —1, is defined through 

(Q\j^0)\p(P,e))=f p m p 8 lM . (3.7) 

Decay p° — > e + e~ allows to measure f p o = 150 MeV, then from isospin 
symmetry we can deduce f p ± = y/2 f p0 = 210 MeV. 

Interesting question: Why do we put T-product of currents in the def- 
inition ()2.1I) of correlators? To answer it we consider the vector current 
correlator H pu in more detail. Lorentz invariance and vector current con- 
servation dictate 

IW<Z)=« Jd i xe i ^(0\T[J^(x)J 1/ (0)}\0)=[q p q I/ -g pu q 2 ] 11(g) . (3.8) 
Inserting 1 in between currents, we obtain 

n(«) = ^5(p-0%o)|(o|J M (o)|x(p)) 2 

9 X(p) 



x J dt p(9°-Po)t + e -»(90+P0)tj _ (3 9) 

From Sokhotsky identity we have 

r-oo i 

/ <ft e ±iat = tt 5(a) ±iV-, (3.10) 
Jo « 
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that give us 
Iml% 2 ) 



-7T ■ 



3g 2 



^5(p-<?)5(p -\q \) (0|J M (0)|X(p)) (3.11) 



As a result we have 



-Iml% 2 

7T 



(3.12) 



with 



p(q 2 )0( qo ) 



"(27T) 



3?2 



£<5< 4 >(g-p)0(po) <0|J M (0)|X(p)} . (3.13) 



Lorentz invariance dictates that p(q 2 ) > 0. Indeed, 

(01^)1X0?)) = ^ + ^^] e"^ 



with p ■ e = 0, and therefore e ■ e = - 
A = 0,1 e. (B = f x m x , see ((577)1) 



(3.14) 

1. From current conservation it follows 



(0|j"(s)|X(p)>(X(p)| Jt(x)|0> = l/xl'm^e 2 = -|/x| 2 ™x < 0. (3.15) 
This gives us 

P(s) = l]|/*| 2 *(«-m^)>0. (3.16) 

Now we can say why we put T-product in correlators - then spectral densi- 
ties, defined only by real particles, are Lorentz invariant and depend only on 
q 2 ! Indeed, we see that the structure e l ^ go ~ po ^ t + e~ t ^ qo+po ^ t in (|3.9|) appears 
due to the presence of the T-product in the definition of the correlator (|3.8|) . 
And just this structure generates as a result 0(|<zo|) m 1(3.120 . 

Relation to the cross section od e+e — — ► hadrons. We have ((3.13)1 for the 
spectral density p{q 2 ). This function naturally appears in 1-photon QED 
description of the process e + e~ — > hadrons if current is the electromag- 
netic current due to quarks: 



ie 

X(p) & u(k)lMk) — (X(p)\Mq)\0) 
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Here k and k! are the momenta of ingoing e + and e . Then we can deduce 

47T0 2 



167r 3 a 2 . . 

^e+e-^hadronsl 5 ) = P\ s ) 



3 s 



R(s), 



(3.17) 



where we explicitly extracted as a factor the cross-section cr e + e - _,.„+„- (s) 
Air a 2 /(3 s) of the process e + e~ — > /x + \f . Equivalently: 



rj/ \ &e+e~— >hadronsV^/ / \ 

R(s) = — ; p(s) 



12 vr 2 



R(s) 



(3.18) 



Now we can look to the quark-hadron duality, that is duality of hadron 
spectral density phad(s), which is measured in r-decay to ^v+hadrons, and 
quark spectral density p per t(s) predicted by QCD (in this case J M is given 
by Eq. (|3.6|) and in the leading order R(s) = N c ): 



S2 f-S2 

PpeA s )ds = / Phad(s)«is 
l J si 



(3.19) 



Looking in Fig. Owe can produce the following observations: 

• Real hadron spectral density is more smooth than in QHO case; 

• Duality is working! 

• Asymptotics starts at s > 3 GeV 2 





\ 

p(s), vector isovector 




1 ALEPH 


J 


\ Ppert(s) ] 

\ \ jjiilUJiU!: 


J 


s [GeV 2 ] 




Fig. 5. We show here perturbative spectral density 47r 2 p pcr t(s) (blue dashed line) in 
comparison with spectral density 4-7r 2 phad(s) measured by the ALEPH Collaboration. 
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4. Condensates in QCD 

In quantum mechanics we saw that in the presence of confinement po- 
tential 



M (r- 1 ) - Mo (r- 1 ) = £ 



6 360 



This difference vanishes at short distances t<1/w and one can calculate 
exact M(fi) perturbatively, expanding in powers of the oscillator potential. 
In QCD confining potential V coni (r) is not even known. How to proceed 
further? The suggestion of QCD SR approach is: 

• To construct perturbative expansion in terms of quark and gluon prop- 
agators; 

• To postulate that quark and gluon propagators are modified by the 
long-range confinement potential; 

• To suppose that this modification is soft: at r — ► the difference 
between exact and perturbative propagators vanishes. 

In realizing this program we write the exact propagator 5 exact (:E, 0) of a 
field ^ as a vacuum average in the exact vacuum |0) 

iS cxact (0,x) = (0|T(V5(x)V>(0))|0). (4.1) 

Wick theorem allows us to write T-product as the sum 



T(V>(x)V>(0)) = if)(x)il>(0)+ :if){x) ^(0) : (4.2) 
of the "pairing" and the "normal" product. Then 

5 exact (0> x) = Q) + (Q| . ^ (x) ^ (Q) . |Q) > (43) 

which can be considered as the starting point to calculate power corrections 
in QCD. The examples in QCD are: (qq) = (0| :q(0) q(0) : |0) referred to 
as quark condensate; (qD 2 q), characterizing average virtuality of the vac- 
uum quarks; gluon condensate (GG) = (Q| (0) (0) : |0), etc. Here 
D[+ = 3^ — igA^ is the covariant derivative and G^ = {i/g)[D^, D u \ is the 
gluonic field strength. 

Condensates and PC AC for pions in QCD. We derived the relations (|3.3|) 
and (|3.4|) ~ ()3.5j) . In order to see their consequences consider now correlator 

IW<7) = i Jd*xe iqx (0\T [j M5 (x)4(0)] |0) = ^Hap^ 2 ) (4.4) 
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and its contraction with q^ 

iq 2 n AP (q 2 ) = -Jd 3 xe-^(0\ [j 05 (0, x); 4(0)] |0> 

-(m u + m d ) Jd i xe igx {0\T[j 5 (x)4(0)] |0> 
= % (uu + dd) +i (m u + m d ) n 55 (g 2 ) . 
Inserting pions in between currents of Hap((? 2 ) in ()4.4|) we have 

fnml U _ ~fn m l 1 



n A p(<r 



m u + m d m 2 — q 2 



q 2^^ mu + md q i 



1 + 

q z 



(4.5) 



(4.6) 



Comparing asymptotics 0(1/ q 2 ) of (|4.5|) and (|4.6|) gives us the famous 
PC AC relation: 



f w m w = - (uu + dd) (m u + m d ) + 0(m q ) . 
In fact we should add other possible PS-meson states to obtain 

fl ml + f 2 , m 2 , + ... = - (uu + dd) (m u + m d ) + Q(m 2 ) 



(4.7) 



(4.8) 



In the chiral limit, m q — > 0, PCAC relation tells us: 

• fn 7^ 0, then — > pion is Goldstone boson; 

• m^i ^ 0, then /„•/ — > =^ no decays 7r ; — ► /Ui^ ! 

• / w ra 130 MeV =^ (gg) « -(260 MeV) 3 at m M = m d = 4 MeV. 

5. QCD SRs for 7r-mesons 

We study axial-axial correlator n^.^g) 

% Jd 4 xe iqx (0\T [j M5 (x)4 5 (0)] |0) = g^Tl^q 2 ) + q^U 2 (q 2 ) 
Hadronic contribution to Borel transform of n^g 2 ): 



<I> 



hadr 



M 2 ) = B Q ^ AP U h 2 ^(q 2 ) 



M 2 M 2 



The following diagrams contributes to the OPE of this correlator 
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$ope (M 2 ) 



$ pert (M 2 ) + $ 2V (M 2 ) 




+ $ 3L (M 2 ) + $ 4Q (M 2 



+ 



$GG (M 2 ) 



with ${2V,3L,4Q}(^ 2 ) = {16,16,144}(vra s (gg) 2 )/(81M 6 ). We see that in 
quark condensate contribution the most important one is <3?4Q. As a result 
we have the following SR for the pion decay constant 



1 + 



7T 



1 (a s GG) 176 Tra s (qq) 2 
12 7r M 2 81 M 4 1 J 



Numerically, as can be seen from Fig. H3 we obtain f n = 0.128 ± 0.13 GeV 
from this SR, whereas in the spectral model with ^4i-meson we obtain 
slightly higher value f w = 0.137 ± 0.13 GeV, to be compared with /^ xp = 
0.132 GeV. 



.136 
0.134 
0.132 
0.13 
0.128 
. 126 
0.124 
0.122 



/^(M 2 ) [ir + continuum] 
Exp. 



sy - 0.69, 



so = 0.63 



s n = 0.57 



» r 2~ 



,1/ 



/„-(M 2 ) [7r + Ai + continuum] 

s = 2.4 _ 



. 135 ■ i r - 



= 2.2 



Exp". "~ — - _ 



AI 2 



3.7 0.8 0.9 



1.1 1.2 1.3 



1.2 1.4 1 . E 



Fig. 6. Left: panel We show here the 1. h. s. of Eq. 15. II . /|(M 2 ), as a function of 
M 2 (red solid line). Right panel: The same as in the left panel, but for the spectral 
model 7r + Ai + continuum. Blue dashed lines correspond to 10% variations of the best 
tresholds. Positions of Ml and My are shown by green vertical arrows in both panels. 
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6. Generalized QCD SRs for mesonic distribution amplitudes 

The pion distribution amplitude (DA) parameterizes the matrix element 
of the nonlocal axial current on the light cone 



(0\d(z) W E(z,0)u(0)\ir(P)) 



z 2 =0 



iUP^Jdx e »( 2P Vr~ W 2 ) • (6-1) 



Here the gauge-invariance is guarantied by the Fock-Schwinger string 



E(z,0) = Vexp 



ig / A^dn* 



in between separated quark fields. The physical meaning of this DA — the 
amplitude of the transition 7r(P) — ► u(Px) + d(P(l — x)). It is convenient 
to represent the pion DA: 



Mx; z" 2 ) = ^0*0 1 + Yl °2n(^ 2 ) Cl / n \2x - 1) 



(6.2) 



n>l 



where C n (2x — 1) are the Gegenbauer polynomials (1-loop eigenfunctions 
of ER-BL kernel) and ip As (x) = 6x (1 — x). This representations means that 
all scale dependence in ^(x) fi 2 ) is transformed to the scale dependence of 
the set |a2(/i 2 ), a^/j, 2 ), . . .}. ER-BL solution at the 2-loop level is also 
possible with using the same representation Q6.2JI _6| |S1 E] • 

In order to construct reliable SRs for the pion DA moments one needs, 
as has been shown in |1U1 lllj . to take into account the nonlocality of QCD 
vacuum condensates. For an illustration of the nonlocal condensate (NLC) 
model we use here the minimal Gaussian model 

{q(o)q(z)) = (qq)e~ lz2lx * /8 . (6.3) 

The single scale parameter X 2 = (k 2 ) characterizes the average momentum 
of quarks in the QCD vacuum and has been estimated in QCD SR approach 
and also on the lattice QH H3 El EH EH : 

\ 2 = 0.35 - 0.55 GeV 2 . (6.4) 

That means that X 2 is of an order of the typical hadronic scale m 2 ~ 0.6 
GeV 2 . We write down as an example the NLC QCD SR for the pion DA 
itself. To produce it one starts with a correlator of currents J^{x) and 
^V5-iv(0) = d(0) (nV) N u(0) with light-like vector n, n 2 = 0, obtains 
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Fig. 7. Left panel: We show here contributions to the R. H. S. of Eq. 16. 5H due to 
perturbative loop (blue dotted line) and due to 4Q-condensate: <p^q(x) - in standard 
QCD SRs, and ^ c (x,M 2 ) (with M 2 = 0.55 - 0.80 GeV 2 ) - in NLC QCD SRs. Right 
panel: Moments (£ )t obtained using the NLC SR 16. 5L are shown by triangles with 
error-bars. We show also for comparison the asymptotic DA moments (squares). 

SRs for the moments (x N ) 7r and then realizes inverse Mellin transform from 
the moments (x N ) n to the DA ip w (x): 



*7T(> 



+ 



8ira s (qq)' 
81M 4 



]T <Pi(x;A). (6.5) 

i=2V,3L,4Q 



The local limit X 2 /M 2 = A —> of this SR is specified by the appearance 
of 5- functions concentrated at the end-points x = and x = 1, for example, 
<Pgg(x; A) = [5(x) + 5(1 - x)] and ip 4Q (x;A) = 9[5(x) + 5(1 - x)}. The 
minimal Gaussian model ()6.3(1 generates the contribution <^4q(x; A) shown 
in the left panel of Fig.[7|in comparison with perturbative one for the stan- 
dard (local) and the NLC types of the SR. We see that in the local version 
due to completely different behaviour of perturbative and condensate terms 
it is difficult to obtain some kind of consistency. Alternatively, the NLC 
contribution is much more similar to the perturbative one — and for this 
reason in the NLC SR we have a very good stability! After processing the 
SR (E3J) for the moments (at fi 2 ~ 1.35 GeV 2 ) 



<Pk(%) (2x — 1) N dx 



(6.6) 



we restore the pion DA cp n (x) by demanding that it should reproduce these 
5 moments and applying the minimally possible number of the Gegenbauer 
harmonics in representation (|6.2j) . It appears the NLC SRs for the pion DA 
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Fig. 8. Left panel: The allowed values of parameters a? and 04 of the bunches 116.711 at 
fi 2 — 1.35 GeV 2 for three values of the nonlocality parameter A 2 : 0.4 GeV 2 (green region), 
0.5 GeV 2 (blue region), and 0.6 GeV 2 (red region). Right panel: The inverse moment 
(x~ 1 ) n , obtained using the NLC SR l|6.5|l . is shown by the green solid line (central value) 
with error-bars, shown as green dashed lines. 

produce a bunch of self- consistent 2-parameter models at fi 2 ~ 1.35 GeV 2 : 



(6.7) 



n=l,2 



,BMS 



-0.130 in the 



The central point corresponds to af MS = +0.188, 
case A 2 = 0.4 GeV 2 , whereas other allowed values of parameters a 2 and 04 
are shown in the left panel of Fig.|H]as the green slanted rectangle ^7j- We 
verify that this solution is self-consistent by estimating the inverse moment 
of the pion DA, (x -1 )^, in two ways. The first is based on Q6.7JI and gives 
us 



(x 



- 1 \ bunch 



3.17 ±0.20. 



The second way uses the special SR for this moment, obtained through the 
basic SR (|6.5|) . It is worth to emphasize here that the moment (x _1 }^ R 
could be determined only in NLC SRs because the end-point singularities 
absent. This SR produces the estimate, see Fig. right panel: 



(x 



-1\SR 



3.30 ± 0.30 



(6.9) 



at ji 2 ~ 1.35 GeV 2 . We see that both estimates are in a good agreement. 

Comparing the obtained pion DA with the Chernyak&Zhitnitsky (CZ) 
one ^H], see Fig. El reveals that although both DAs are two-humped they 
are quite different: our DA is strongly end-point suppressed. This can also 
be verified in the right panel of the figure, where contributions of different 
bins to inverse moment (x -1 ),,-, calculated as f x+a02 (fi(x)dx and normalized 
to 100%, are shown for CZ and BMS DAs. 
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0.4 O.G 

X 

Fig. 9. Left panel: We show here the comparison of curves for three DAs — BMS 
(green solid line), CZ (red dashed line), and the asymptotic DA (black dotted line). 
Right panel: Histograms for contributions of different bins to inverse moment (x^ 1 )^ 
are shown for CZ and BMS DAs. 



7. LCSR analysis of CLEO data on E y ~ f * 7V (Q 2 ) and pion DA 

Why does one need to use Light-Cone SRs (LCSRs) in analyzing the ex- 
perimental data on 7* (Q)^(q) — > 7r°-transition form factor? For Q 2 S> m 2 , 
q 2 <C m 2 the QCD factorization is valid only in the leading twist approxi- 
mation and the higher twists are of 
importance [19] . The reason is evi- 
dent: if — ^ one needs to take 
into account interaction of a real 
photon at long distances of order 
of 0(1/ y/q 2 ). To account for long- 
distance effects in perturbative QCD 
one needs to introduce the light-cone DA of the real photon. 

Instead of doing so, Khodjamirian [2U] suggested to use the LCSR ap- 
proach, which effectively accounts for long-distances effects of the real pho- 
ton using the quark-hadron duality in the vector channel and dispersion 
relation in q 2 . 

We refined the NLO analysis of the CLEO data [U by taking into ac- 
count the following items: 

(i) an accurate NLO evolution for both tp(x, Q 2 '') and a s (Q 2 ) with ac- 
counting for quark thresholds; 

(ii) the relation between the "nonlocality" scale and the twist-4 magnitude 

~ ^q/ 2 was useci to re-estimate £f w _ 4 = 0.19 ±0.02 at X 2 = 0.4 GeV 2 ; 

(iii) the possibility to extract constraints on (x^ 1 )^ from the CLEO data 
and to compare them with what we have from NLC QCD SRs. 

The results of our analysis [22] are displayed in Fig. 1101 Solid lines in 
all figures enclose the 2<r-contours, whereas the ler-contours are enclosed 
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Fig. 10. Three 2(j-contours of the admissible regions following from the analysis of 
the CLEO data for different values of S 2 : (a) for A 2 = 0.6 GeV 2 and <5t w _ 4 ~ 
(0.29 ± 0.03) GeV 2 ; (b) for A 2 , = 0.5 GeV 2 and 4 W _ 4 = (0.235 ± 0.025) GeV 2 ; (c) for 
A 2 = 0.4 GeV 2 and = (0.19 ± 0.02) GeV 2 . 



by dashed lines. The three slanted and shaded rectangles represent the 
constraints on (02, 04) posed by the QCD SRs \Tj\ for corresponding values 
of A 2 = 0.4, 0.5, 0.6 GeV 2 (from left to right). All values are evaluated at 
/x 2 = 2.4 GeV 2 after the NLO evolution. 

We see that the CLEO data definitely prefer the value of the QCD 
nonlocality parameter A 2 = 0.4 GeV 2 . We also see in Fig. HUT c) (and this 
conclusion was confirmed even with 20% uncertainty in twist-4 magnitude, 
see also Fig. Ill[) that CZ DA (■) is excluded at least at 4<r-level, whereas 
the asymptotic DA (♦) — at 3<7-level. In the same time our DA (JC) and 
most of the bunch (the slanted green-shaded rectangle around the symbol 
X) are inside lcr-domain. Instanton-based models are all near 3u-boundary 
and only the Krakow model |23| . denoted in Fig. ^2 by symbol ♦, is close 
to 2<r-boundary. 

In the left panel of Fig. ^2 we demonstrate the la-, 2a- and 3cr-contours 
(solid, dotted and dashed contours around the best-fit point (+)), which 
have been obtained for values of the twist-4 scale parameter <^ w _ 4 = [0.15 — 




Fig. 11. The results of the CLEO data analysis for the pion DA parameters ((a; -1 )^ xp /3 — 
1, evaluated at /z 2 , fa 1 GeV 2 , in the left panel, and a 2 and 04, evaluated at — 
5.76 GeV 2 - in the right panel). In the right panel the lattice results of |29| are 
shown for comparison as shaded area, whereas the renormalon-based la-ellipse of |27| is 
displayed by the green dashed line. 
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0.23] GeV 2 . As one sees from the blue dashed line within the hatched band, 
corresponding in this figure to the mean value of (x _1 )| R /3 - 1 and its 
error bars, the nonlocal QCD sum-rules result with its error bars appears 
to be in good agreement with the CLEO-constraints on (x -1 )^ xp at the ler- 
level. Moreover, the estimate (x _1 )| R is close to (x _1 )^ M /3 - 1 = 0.24 ± 
0.16, obtained in the data analysis of the electromagnetic pion form factor 
within the framework of a different LCSR method in |241 125 j . These three 
independent estimates are in good agreement to each other, giving firm 
support that the CLEO data processing, on one hand, and the theoretical 
calculations, on the other, are mutually consistent. 

Another possibility, suggested in |26| . to obtain constraints on the pion 
DA in the LCSR analysis of the CLEO data - to use for the twist-4 con- 
tribution renormalon-based model, relating it then to parameters ai and 
04 of the pion DA. Using this method we obtain [27] the renormalon-based 
constraints for the parameters a>2 and 04, shown in the right panel of Fig. 
1111 in a form of lcr-ellipses (dashed contour). 

New high-precision lattice measurements of the the pion DA second 
moment (£ 2 )tt = Jq(2x — l) 2 fTr(x) dx appeared rather recently [2*51 12*""]. 
Both groups extracted from their respective simulations, values of 02 at the 
Schmedding-Yakovlev scale /i| Y around 0.24, but with different error bars. 

It is remarkable that these lattice results are in striking agreement with 
the estimates of ci2 both from NLC QCD SRs 17 and also from the CLEO- 
data analyses — based on LCSR — |21M22j . as illustrated in the right panel of 
Fig. 111! where the lattice results of [!"""j] are shown in the form of a vertical 
strip, containing the central value with associated errors, being smaller than 
in [2H] . Noteworthily, the value of 02 of the displayed lattice measurements 
(middle line of the strip) is very close to the CLEO best fit in [21] (+). 

8. Pion form factor and CEBAF data 

It is worth to mention here the results of our analysis of the pion electro- 
magnetic form factor using NLC dictated pion DA and Analytic Perturba- 
tive QCD "SUJ. These results are in excellent agreement with CEBAF data 
on pion form factor, shown as diamonds in the Fig. 1121 where the green strip 
includes both the NLC QCD SRs uncertainties, generated by our bunch of 
the allowed pion DA, and by the scale-setting ambiguities at the NLO level. 

From the phenomenological point of view, the most interesting result 
here is that the BMS pion DA [T7j (out of a "bunch" of similar doubly- 
peaked endpoint-suppressed pion DAs) yields to predictions for the elec- 
tromagnetic form factor very close to those obtained with the asymptotic 
pion DA. Conversely, we see that a small deviation of the prediction for 
the pion form factor from that obtained with the asymptotic pion DA does 
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Fig. 12. Predictions for the scaled pion form factor calculated with the BMS bunch (green 
strip) encompassing nonperturbative uncertainties from nonlocal QCD sum rules |17| and 
renormalization scheme and scale ambiguities at the level of the NLO accuracy. The 
dashed lines inside the strip indicate the corresponding area of predictions obtained with 
the asymptotic pion DA. The experimental data are taken from 31 (diamonds) and \'S2\. 
(triangles). 

not necessarily imply that the underlying pion DA has to be close to the 
asymptotic profile. Much more important is the behavior of the pion DA in 
the endpoint region x — * , 1. 



Let me conclude with the following observations: 

• NLC QCD SR method for the pion DA gives us the admissible bunches 
of DAs for each value of \ q . 

• NLO LCSR method produces new constraints on the pion DA param- 
eters (02 and 04) in conjunction with the CLEO data. 

• Comparing results of the NLC SRs with new CLEO constraints allows 
to fix the value of QCD vacuum nonlocality: A 2 - ~ 0.4 GeV 2 . 

• This bunch of pion DAs agrees well with recent lattice data and with 
JLab F(pi) data on the pion form factor. 

I also suggest to the reader to look in the very interesting discussion of the 
QCD SR approach and its developments written by one of its creators |34| . 



9. Conclusions 
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